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We calculated the energy spectra and the fluxes of electron neutrinos emitted in the process 
of the evaporation of primordial black holes (PBHs) in the early universe. It was assumed that 
PBHs are formed by a blue power-law spectrum of primordial density fluctuations. We obtained 
the bounds on the spectral index of density fluctuations assuming validity of the standard picture of 
gravitational collapse and using the available data of several experiments with atmospheric and solar 
neutrinos. The comparison of our results with the previous constraints (which had been obtained 
using diffuse photon background data) shows that such bounds are quite sensitive to an assumed 
form of the initial PBH mass function. 
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Some recent inflation models (e.g., the hybrid inflationary scenario p|) predict the "blue" power-spectrum of 
primordial density fluctuations. In turn, as is well known, the significant abundance of primordial black holes (PBHs) 
is possible just in the case when the density fluctuations have an n > 1 spectrum [n is the spectral index of the initial 
density fluctuations, n > 1 spectrum is, by definition, the "blue perturbation spectrum"). 

Particle emission from PBHs due to the evaporation process predicted by Hawking Q may lead to observable 
^•f-) , effects. Up to now, PBHs have not been detected, so the observations have set limits on the initial PBH abundance 
■ or on characteristics of a spectrum of the primordial density fluctuations. In particular, PBH evaporations contribute 
' to the extragalactic neutrino background. The constraints on an intensity of this background (and, correspondingly, 
on an PBH abundance) can be obtained from the existing experiments with atmospheric and solar neutrinos. The 
obtaining of such constraints is a main task of the present paper. 

The spectrum and the intensity of the evaporated neutrinos depend heavily on the PBH's mass. Therefore, the 
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' great attention should be paid to the calculation of the initial mass spectrum of PBHs. We use in this paper the 
, following assumptions leading to a prediction of the PBH's mass spectrum. 

1. The formation of PBHs begins only after an inflation phase when the universe returns to the ordinary radiation- 
dominated era. The reheating process is such that an equation of state of the universe changes almost instantaneously 
into the radiation type (e.g., due to the parametric resonance J| ) after the inflation. 

2. It is assumed, in accordance with analytic calculations [[|||] that a critical size of the density contrast needed for 
+3 . the PBH formation, 5 C , is about 1/3. Further, it is assumed that all PBHs have mass roughly equal to the horizon 
^ ' mass at a moment of the formation, independently of the perturbation size. 

3. Summation over all epochs of the PBH formation can be done using the Press-Schechter formalism ||. This 
formalism is widely used in the standard hierarchial model of the structure formation for calculations of the mass 

k> , distribution functions (see, e.g., J7|). 

It was shown recently that near the threshold of a black hole formation the gravitational collapse behaves as a 
critical phenomenon [|S|. In this case the initial mass function will be quite different from that follows from the 
standard calculations of refs. ^Q.The first calculations of the PBH initial mass function (e.g., in ref. ||) have been 
done under the assumption that all PBHs form at the same horizon mass (by other words, that all PBHs form at 
the smallest horizon scale immcdiatly after reheating). The initial PBH mass spectrum for the case of the critical 
collapse, based on the Press-Schechter formalism, was obtained in the recent work ||. 

The calculations in the present paper are based on the standard [Qj5| picture of the gravitational collapse leading 
to a PBH formation. The case of the critical collapse will be considered in a separate work. 

The plan of the paper is as follows. 

In Sec.[n] we give, for completeness, the brief derivation of a general formula for the initial PBH mass spectrum. 
The final expression is presented in a form, which is valid for an arbitrary relation between three physical values: 
the initial PBH mass, the fluctuation mass (i.e., the mass of the perturbed region) at a moment of the collapse, 
and the density contrast in the perturbed region (also at a moment of the collapse). This expression contains the 
corresponding results obtained in refs. [^|^] as particular cases. As in refs. the derivation is based on the linear 

perturbation theory and on the assumption that a power spectrum of the primordial fluctuations can be described by 
a power law. 



In Sec. Ill we derive the approximate formula for a calculation of the extragalactic neutrino background from PBH 



evaporations. We do not use cosmological models of the inflation and of the spectrum of primordial fluctuations, so 
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there are two free parameters: the reheating temperature and the spectral index. At the end of the section, some 
examples of instantaneous neutrino spectra from the evaporation of an individual black hole are presented. 

In Sec. IV we determine the explicit dependence of an background intensity on the spectral index (by normalization 
of standard deviation of the density contrast at horizon crossing on COBE data). After this, numerical calculations 
of the neutrino background become possible. For background neutrino energies ~ 10 — 100 MeV we studied relative 
contributions to the background intensity of different cosmological redshifts (and it is shown that at high reheating 
temperatures the characteristic values of the redshifts are very large). Further, several results of numerical calculations 
of the neutrino background spectra are presented. Possibilities of constraining the spectral index using experiments 
with neutrinos of natural origin are discussed. 

In Sec.|v| the spectral index constraints, followed from our calculations and from available data of the neutrino 
experiments, are given. It is proved that, almost in all intervals of the reheating temperatures considered in this work, 
effects of an background neutrino absorption in the space are small and can be neglected. 



II. THE INITIAL MASS SPECTRUM OF PBHS 



As it pointed out in the Introduction, we use the Press-Schechter formalism which allows to carry out the 
summation over all epochs of PBHs formation. According to this formalism, the mass spectrum of density fluctuations 
( i.e., the number density of regions with mass between M and M + dM ) is calculated by the formula 



n (M, S c ) dM = — 



dM y ' c ' 



dM . (2.1) 
Here, (3{M, S c ) is the fraction of regions having sizes larger than R and density contrast larger than S c , 

(2.2) 
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f3(M, S c ) = 2 / P(M, S)d6 



P(M, S) = 
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2tT<Tr(M) 
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2o\ (M) 
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Here, S is the initial density contrast, ur is the standard deviation of the density contrast of the regions with size R 
and mass M . 

It is convenient to introduce the double differential distribution n(M,S), the integral over which gives the total 
number of fluctuated regions, 

rr r 

n(M, S)dMdd , n(M, S c ) = / n(M, S)dd . (2.4) 

Using Eqs.( [2.l| - |2.3| ), for n(M, 6) one has the expression 

It ft 1 da R 



n(M, 6) 



it M a 2 R (M) dM 



a R (M) 



- 1 



exp 



2a\{M) 



(2.5) 



To obtain the mass spectrum of PBHs one must introduce the Mbh variable. Besides, we will use the variable 5' , 
connected with 5 by the relation 



5' 



M 
Ml 



2/3 



(2.6) 



Here, Mj is the horizon mass at the moment of a beginning of the growth of density fluctuations. This new variable 
is the density contrast at the moment of the collapse. The new distribution function is 



,„ r r/ . r , d5 dM 

nBH\M B H,0 ) = n(M B H,0) — - 



dS' dM BH 

Further, we assume that there is some functional connection between Mbh, M and S': 

M BH = f(M,5') . 



(2.7) 



(2.8) 
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In this case one can rewrite Eq.(2.7) in the form: 

tibh{M bh , 6') = n | M, 6' (j£ 



-2/3 N 



M 

M 



-2/3 



1 



df(M, 5')/dM 



and the PBH mass spectrum is given by the integral 

n B H{M BH ) = J n BH {M BHl 8')d8' . 
Now, to connect the mass spectrum with the spectral index one can use the relations 



or = o H (M) 



M 



-2/3 



o H {M) ~ M-ir 



(see Sec. IV for details). Using Eqs.(2.11) in the expression (|2.5|) for n(M,5) one obtains: 



n(M, 8) 



n + 'i 2 pi 1 



M 



-4/3 



Substituting this expression in the r.h.s. of the Eq.(B7|) and using Eq(2.10) one has: 



ti B h{M B h) = 



n + 3 



M 2 



oh 



8' 2 



- 1 



1 



df(M, S')/dM 



dS' 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



Values of M in r .h.s of Eq.( 2.13 ) are expressed through M BB , 8' by the relation (|2 

Formula ( 2.1 3| ) is the final expression for the PBH mass spectrum, the main result of this Section. It is valid for 
any relation between PBH mass M BB , mass of the original overdense region M and the density contrast 8' . 

In the particular case of a near critical collapse 



and 



f(M, 8') = M- /3 M 2/3 k(8' - J c ) 7fc = ^mV 3 M 2//3 



M = M^M" 1/2 r 3/2 ; df/dM = W^M^e 12 



(2.14) 



(2.15) 



Eq.(2.14) can be rewritten in the form, derived in the ref. p[, using the connection between M and the horizon 
mass Mh (which is equal to a fluctuation mass at the moment when the fluctuation crosses horizon): 



M h = M-'°M 



l/3 M 2/3 



The resulting formula for the spectrum 
n{M BH ) 
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S c OH 



V2 



' H 
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e ^£ 3/2 dS' , 



(2.16) 



(2.17) 



coincides, as can be easily proved, with the expression derived in ref. ||. 

The mass spectrum formula for the Carr-Hawking collapse can be obtained analogously, with the relation f(M, 8') 

7 1 / 2 M i L/3 Af 2 / 3 , or from Eq.( |2.17| ) using the substitutions 

„V2 



and the approximate relation 



1 ( 8' 2 

T H 
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8 C -> 7 



7e 



In this case one obtains the expression obtained in ref. 

n + 3 Fl 



n BH (M BH ) = 



i 7/i pM ,2 M B T^ 



exp 



2a% 



(2.18) 
(2.19) 

(2.20) 



One can see from Eqs.(§T7|) and (|2~2p| ) that the PBH mass spectrum following from the Press-Shechter formalism 

/ -5/' ' 

has quasi power form in both considered cases l~ M BB 



In Carr-Hawking case M 



in in 
BH 



Mh , in contrast with 



this in the critical collapse case it is possible that M B B <SC Mh ■ However , the low mass part of the PBH spectrum 
is suppressed by the factor £ 3 / 2 in the integral in Eq. (2.17). 
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III. NEUTRINO DIFFUSE BACKGROUND FROM PBHS 



The starting formula for a calculation of the cosmological background from PBH evaporations is |lC 
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S(E) = I n comT -^f (E(l + z)) dV com . (3.1) 



Here , n com is the comoving number density of the sources (in our case the source is an evaporating PBH of the 
definite mass m ), ao is the scale factor at present time, t — to, f(E) is a differential energy spectrum of the source 
radiation, V com is a comoving volume of the space filled by sources, therefore 

dV C o m = al— f dp dVt . (3.2) 

v 1 - k P 

Here, k is the curvature coefficient, and p is the radial comoving coordinate. Using the change of the variable, 

dp dt 



y/1 - kp 2 a ' 

the comoving number density can be expressed via the initial density rii, 



(3.3) 



p,, s (M^(^) 3 (£) = »*(^) ■ (3-4) 



(3.5) 



Substituting Eqs. (3.2) - (3.4) in Eq. (3.1) one obtains 

S(E)= ni f dt^(^) f(E(l + z)) . 
J a \a J 

In our concrete case the source of the radiation is a Hawking evaporation: 

n l f(E(l + z))=[dmn BH (m,t)f H {E(l + z),m) . (3.6) 



Here , riB^im, t) is the PBH mass spectrum at any moment of time, fn(E, m) is the Hawking function Eq], 

f tw \ 1 r s (E,m) 

/g( ^ TO)= 2^ e , P (^)- M p ' (3J) 

Here, T s (E,m) is the coefficient of the absorption by a black hole of a mass m, for an particle having spin s and 
energy E. 

Initial spectrum of PBHs is given by Eq. (2.20). The minimum value of PBH mass is 7 M i , so we must add to 

1 /2 

the initial spectrum expression the step factor Q{ni BH — 7 MJ . The connection of the initial mass value Mbh 
and the value at any moment t is determined by the solution of the equation |Tl| ] : 

dm = a(m) 
dt m 2 

The function a(m) accounts for the degrees of freedom of evaporated particles and determines the lifetime of a black 
hole. In the approximation a = const the solution of Eq.( |3.8| ) is: 

M BH = {Sat + to 3 ) 1/3 . (3.9) 

This decrease of PBH mass leads to the corresponding evolution of a form of the PBH mass spectrum. At any moment 
one has 



2 r 

n BH (m, t)dm = ^n BH ((3ai + to 3 ) 1 / 3 ) x 6 to - (( 7 1 / 2 M l ) 3 

(Sat + to 3 ) ' ^ ' L v 



Zat 



1/3" 



dm. (3.10) 
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Substituting Eqs. (3.7), (3.10) in the integral in Eq.(3.6) , we obtain the final expression for the spectrum of the 
background radiation: 



Art J a \a 



dm 



{Zat + m 3 ) 1 / 3 







/2 m,: 



3at 



nsH 



1/3' 



[Bat + m c 



,!/3 



(3.11) 



f H (E(l+z),m) 



One should note that the corresponding expressions for the spectrum in refs. jl2],[L3) contain the factor (^-) 3 instead 

of the correct factor ffj-) ■ It leads to a strong overestimation of large z contributions in S(E) (see below, Fig. 3). 
It is convenient to use in Eq. (3.11) the variable z instead oft. In our case = SIk = 0) we have 



dt 

dz 



{n mo (i + z) 3 + n ro (i + z) 4 ) 1/2 



H (l + z) 

n ro = (2.25 • lo^h 2 )- 1 , n mo = i - n ro 

The factor (f 1 ) 3 can be expressed through the value of t eq : 



(3.12) 



(i + z eq y 



eq 



3/2 



H- 3/2 (2.25 ■W 4 h 2 y m t 3 i /2 



(3.13) 



Integrating over PBH's mass in Eq. (3.11) one obtains finally, after the change of the variable ( on z, the integral 
over z: 



S(E) = J dlog 10 (z + l)F(E,z) 



(3.14) 



In analogous calculations of the photon diffuse background integral over z in the expression for S(E) is cut off at 
z = zo ~ 700 because for larger z the photon optical depth will be larger than unity H]. In contrast with this , 
interactions of neutrinos with the matter can be neglected up to very high values of z. Therefore the neutrino diffuse 
background from PBH evaporations is much more abundant. The neutrino absorption effects are estimated below, in 
Sec.0 

The evaporation process of a black hole with not too small initial mass is almost an explosion. So, for a calculation 
of spectra of evaporated particles with acceptable accuracy it is enough to know the value of a for an initial value of 
the PBH mass only. Taking this into account and having in mind the steepness of the PBH mass spectrum, we use 
the approximation 



a(m) = a(M] 



nun \ 
BH ) 



a(-f 1/2 Mi) , 



(3.15) 



and just this value of a is meant in the expressions (|3.9D -( |3.1l|) . 

The very detailed calculation of the function a(m) was carried out in the works |l5|,[l6|]. Here we use the simplified 
approach in which a{m) is represented by the dependence 



a = a a + ai ■ 10 25 • <d(bi - logio(m)). 



(3.16) 



TABLE I. Coefficients a t and b t in Eq.(|3.16|) 



i 




u,d 
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r 


c 


b 


W ± 


Z° 


t 


(Li 


3.12 


9.36 


9.36 


4.96 


3.12 


9.36 


9.36 


1.86 


0.93 


9.36 


U 


14.47 


14.01 


13.77 


13.67 


13.23 


13.20 


12.79 


11.57 


11.55 


11.27 



Here, &(x) is the Heaviside step function. Coefficient ao gives the summary contribution of e + , e , v and 7 and is 
equal to 8.42 • 10 25 g 3 sec~ 1 All other coefficients are collected in the Table @. 
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The coefficients bi determine the value of the PBH mass beginning from which particles of i-type can be evaporated. 



This value, M BH , is obtained from the relations 
So, 

bi = log w M b BH . (3.18) 
The resulting function a(m) is shown on Fig.^ For comparison, the corresponding dependence from the work |l(f | 



(drawn using th e Eq. 7 of |16|) is also shown. 

The formula (3.11), as it stands, takes into account the contribution to the neutrino background solely from a 
direct process of the neutrino evaporation. If we suppose that particles evaporated by a black hole propagate freely, 
the calculation of other contributions to the neutrino background can be performed using our knowledge of particle 
physics pxlq]. 

On Fig!p the typical result of our calculation of instantaneous neutrino spectra from evaporating black hole is 
sho wn. The spectrum of the straightforward (direct) (y e + v e ) - emission is given by the Hawking function Jh{E, to), 
Eq.( |3.7[ ). The (y e + v e ) spectrum arising from decays of (/x + + /i~) evaporated directly is calculated by the formula 

f M (E,m) = J f H (E„ m ) dn ^ E) dE, , (3.19) 

where dn^/dE is the neutrino spectrum in a [i - decay. To evaluate the electron neutrino spectrum resulted from 
fragmentations of evaporated quarks, we used the simplest chain: 

(u, d) — quarks — * ir — > fi — > v e , (3.20) 

and the formula 

f { q)(T? s f f fP , dnUQ d£ dnl(E^E v ) dn^j^E) 

f^(E,m) = J f H{Eq ,m)^ r ■ — ■ -—dE^dE, . (3.21) 

Here, dn\jd^ is the q — > n fragmentation function, for which the simple parametrization was taken: 

and dn^/dEft is the /i - spectrum in a decay ir — > fi + v^. 

Analogous formulas are used for calculations of the neutrino spectrum from other channels of the neutrino produc- 
tion, for instance from the decays of evaporated VF-bosons (W — > e + v e , W — > \i — > v e ). 

The relative contribution of different channels to the total v e spectrum depends on the black hole temperature. 
One can see from Fig. 2, that at high temperature decays of massive particles evaporated by the black hole become 
very important at high energy tail of the spectrum (if the corresponding branching ratios are not too small) . 

The total instantaneous neutrino spectrum from a black hole evaporation is given by the sum 

f(E,m)=f H (E,m) + fM{E,m) + fM(E,m) + ... , (3.23) 
and the total background neutrino spectrum is given by the same Eq.( |3.11 ), except the change fjj(E,m) — > f(E,m). 



IV. CONSTRAINTS ON THE SPECTRAL INDEX 

The spectral index of initial density fluctuations is defined by the relations 

-3_ 1 f |r|3w2, 



141= Ak* 



(4.1) 



G 



FIG. 1. The parameter a in the (3.16) for the rate of the PBH mass loss, as a function of the PBH mass. Solid line is the 
present result, dashed line is drawn using Eq.(7) of Ref. [16]. 
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FIG. 2. Instantaneous electron neutrino spectra from a PBH evaporation, for the several channels of the neutrino production. 
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Here, 6k and Wk are Fourier transforms of the density field 8(x) and the window function of comoving size r, 
respectively, V r is the effective volume filtered by W r . 



One obtains from Eq.(4.1 



i.3 
K fl 



where kfi is the comoving wave number, characterizing the perturbed region, 



kfi = 



R 



(4.2) 



(4.3) 



and i? is the physical size of this region at arbitrary moment t. Using the connection of kfi with the fluctuation mass 



(4.4) 



and introducing the horizon mass M/, (w hich is equal to the fluctuation mass at the moment when the fluctuated 
region crosses horizon) we can rewrite Eq.(4.1) in the form: 



( M hor (t) \ V \ 

Here, Mh or is the horizon mass at t , aniM^) is the standard deviation at horizon crossing. 
At the initial moment of time one has : 

M hor {ti) = Mi , M{ti)=M . 
The form of the cr^(M^) function depends on the spectral index: 



(4.5) 



Off 



Ml 



radiation dominance 
matter dominance . 



Eqs.( Q 1 ) in Sec|n]are obtained from Eqs. (i]|)-(0). 

From COBE data we know the normalization of aij\Mh) at present horizon size [f^Jl9[ |: 



a H {M, 



h,0 



9.5 • 10 



-5 



h ,0 



(4.6) 



(4.7) 



Now one can easily show that crjy(M^) for radiation dominance case is connected with aniMhfi) by the following 
approximate relation: 



(M h ) = a H (M h , ) (J^>° 



eq 



M h , 



eq 



(4.8) 



Our calculation of neutrino spectra from evaporating PBHs contains two parameters: a spectral index n and a time 
of an end of the inflation tj (which, by assumption, is a time when density fluctuations develop). We assume that at 
ti the universe has as a result of the reheating the equilibrium temperature Trh- The connection of Tru and ti is 
given by the standard model (in formulas of this section we use the convention h = c = 1): 



U = 0.301g s 



-1/2 M p i 



0.24 



T 2 

RH 



T 2 RH (MeV) 



(4.9) 



(<7* ~ 100 is the number of the degrees of freedom in the early universe). 

On Figjl the integrand of the background spectrum integral (Eq.(3.14)) is shown as a function of redshift z for 
several values of the parameter Trh- Each curve on Fig.|| has a strong cut off on some redshift value. This feature is 
connected with the existence in our model the minimum value of PBH mass, 



BH =1 ' M i 



(4.10) 







FIG. 3. Redshift dependence of the of the integrand in the expression ( 3.14 ) for a neutrino background spectrum, for two 
values of the neutrino energy and for several values of the parameter Trh ■ 
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The PBH mass spectrum is steeply falling function of the mass, so the masses near minimum give a largest contribution 
to the neutrino background. The moment of their evaporation is, approximately, 

tev „ ^MUL , (4.11) 
6a 

and the corresponding redshift is determined by the relation 

ft \ 1/2 

z ev + 1 » (z eq + 1) ( -2* ) . (4.12) 

Larger masses evaporate at larger times and smaller redshifts. If, for instance, Trh = 10 10 GeV, one has ti = 
0.24 • 10~ 26 s ; Mi = 7, 5 • 10 10 g , z ev ~ 10 7 . So, at Trh = 10 10 GeV the redshifts of order of 10 7 give a largest 
contribution to the neutrino background, and this is clearly seen at Fig^. 
The cut off value z ev strongly depends on Trh'- 

zev~~^=~ (M™")" 3/2 ~ t; 3/2 ~ T RH . (4.13) 

V 

On Figs.^ spectra of electron neutrino background are shown separately for several channels of the neutrino pro- 
duction. One can see from these figures that the contribution to the summary background from the direct v e emission 
is dominant (at least, for E > 10 MeV) and the relative importance of different channels changes with an increase of 
Trh- It is seen also that the contribution of the quark fragmentation channel is very small at E > 10 MeV (< 1%), 
so the evident underestimation of this channel in our calculation, connected with the neglect of the contribution of 
heavy quark fragmentations, has no particular importance. 

Some typical results of electron neutrino background calculations are shown on Figs.^,^ The main features of such 
spectra have been revealed in previous works | |12|,p6|p7| | : E~ 3 dependence above 100 MeV, the flattening out of the 
spectra at lower energies and the shift of the turnover energy with a change of the Trh value p2| . According to 
Figs.|||3]the turnover energy is about 1 MeV at Trh = 10 9 GeV. 

For an obtaining of the constraints on the spectral index we use three types of neutrino experiments. 

1. Radiochemical experiments for the detection of solar neutrinos. There are data from the famous Davis experiment 
p0| and the Ga — Ge experiment [ gjj . The cross section for the neutrino absorption via a bound-bound transition 
was calculated using the approximate formula 

^((l) 2 +(2r) 2 (<?) 2 )peE e , E<l00MeV 



(J=< , n i y gv j \«/ j^e , ( 4M ) 

const. , E > lOOMeV . 

In the case of the Davis experiment (CI — Ar reaction) we take into account the super-allowed transition only, for 
which 

<1) 2 = 3 , (a) 2 = 0.2 , E thr » 5MeV. 

In the Ga — Ge case the main contribution gives the ground state - ground state transition (Ethr = 0.242MeV). The 
corresponding cross section is, for E < 100MeV\ 

a (v + Ga -► Ge + e) = 0.646 • 10~ 44 p e £ e - (4.15) 

The number of target atoms N T is about 2.2 • 10 30 for the CI — Ar experiment |EQ] and ~ 10 29 for the Ga — Ge 
experiments pj}| . The average statistics is ~ 1.5 day^ 1 (CI — Ar) and ~ 1 day^ r ^(Ga — Ge). The constraint is 
calculated using the relation 



4tt- N t 



10 5 • J S(E)a(E)dE < 1. (4.16) 



2. The experiment on a search of an antineutrino flux from the Sun [ p2| . In some theoretical schemes (e.g., in 
the model of a spin - flavor precession in a magnetic field) the Sun can emit rather large flux of antineutrinos. LSD 



experiment |22| sets the upper limit on this flux, ^y/^ u < 1.7%. In this experiment the neutrino detection is carried 



out using the reaction 

v e +p — > n + e + . (4-17) 
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FIG. 4. Separate contributions to neutrino background spectra for different channels of the neutrino production during PBH 
evaporation, for fixed values of the parameters Trh and n. 
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FIG. 5. Total energy spectra of electron neutrino background calculated with n — 1.25 and with different values of the 
parameter Trh- 




FIG. 6. Total energy spectra of electron neutrino background calculated with Trh = 10 9 GeV and with different values of 
the spectral index. The theoretical atmospheric neutrino spectrum at Kamiokande site [] (avereged over all directions) is shown 
by the dashed line. 




The number of target protons is ~ 8.6 • 10 per 1 ton of the scintillation detect or, a nd the obtained upper limit is 



0.28 antineutrino events per year per ton [|22|. The cross section of the reaction ( 4.17 ) is well known (see, e.g., p3[). 
It grows with the neutrino energy up to E ~ 2 GeV , and is a constant (~ 0.5 • 10~ 38 cm 2 ) at larger energies. The 
product of this cross section and a PBH antineutrino background spectrum has the maximum at E ~ lOOMeV. The 
constraint is determined from the condition that the calculated effect in the LSD detector is smaller than the upper 
limit obtained in ref. |22| 

3. The Kamiokande experiment on a detection of atmospheric electron neutrinos pi}] . In this experiment the 
electrons arising in the reaction 

vf m + n -» p + e~ (4.18) 

in the large water Cherenkov detector were detected and, moreover, their energy spectrum was measured. This 
spectrum has a maximum at the energy about SOOMeV. The spectrum of the atmospheric electron neutrinos is 
calculated (see, e.g., [p5| ) with a very large accuracy (assuming an absence of the neutrino oscillations) and the 
experimentally measured electron spectrum coincides, more or less, with the theoretical prediction. The observed 
electron excess at E ~ lOOAleV (which is a possible consequence of the oscillations) is not too large. We use the 
following condition for an obtaining the our constraint : the absolute differential intensity of the PBH neutrino 
background at the neutrino energy E ~ O.SGeV cannot exceed the theoretical differential intensity of the atmospheric 
electron neutrinos at the same energy (otherwise the total electron energy spectrum is strongly different from the 
observed one). 



V. RESULTS AND DISCUSSIONS 

Fig.0 shows our results for the spectral index constraints. It is seen that the best constraints are obtained using 
the Kamiokande atmospheric neutrino data and the LSD upper limit on an antineutrino flux from the Sun. 

The behavior of the constraint curve n(Tftjy) is sharply different from that was obtained in ref. |l3f| (authors of ref. 
p3[ used the initial PBH mass function following from the near critical collapse scenario with a domination of the 
earliest epoch of PBHs formation and diffuse extragalactic photon data). 

The slightly non-monotonic behavior of the curves Figj^ near Trh ~ 3 TO 8 GeV is connected with the sharp increase 
of the function a(m) atmss 10 14 g due to a switching on the quark evaporations (see Fig{j]). Indeed, according to the 
relations Eq.( 4.13| ), the minimum value of PBH mass in the initial mass spectrum is inversely proportional to T^ H . 



More exactly, the relation between MgjSp and Trh is 



7.2 ■ 10 30 
TjjGeV) 



M bh(9) = 7r7i • (5-1) 



Thus, when the param eter Trh passes over the value ~ 3 • 10 s GeV, M^^p becomes smaller than 10 14 g, and the 



value of a used in Eqs.( 3.10 - 3.11 ) increases on ~ 20 units. 

One can compare our spectral index constraints with the corresponding results of ref. |12j , where the same initial 
mass spectrum of PBHs had been used. Some difference in resulting constraints (An ~ 0.02) may be connected simply 
with the fact that authors of ref. [Q use slightly different formula for <7Ef (M/i), namely, 

1 — n 

o-ff(M,,) = a H (M h . ) (-^) 4 ■ (5.2) 

\Mh ,0/ 

In other respects the constraints are quite similar although in ref. |T^ | they were based on diffuse extragalactic photon 
background data. 

Constraints on the spectral index, obtained using the neutrino background calculations, are sensitive to a possible 
existence of a QCD photosphere (2|j2^] around black holes only in one case: if the density of the outward-propagating 
plasma is so high that it is opaque for neutrinos. It was shown above (see fig^) that the fragmentations of evapo- 
rated quarks are inessential and can be neglected when calculating the electron neutrino background (in spite of the 
fact that neutrinos from these fragmentations are very abundant (fig.^J)). The reason is simple: neutrino energy is 
strongly redshifted during cosmological expansion. It is clear that the fragmentations of photosphere's quarks are also 
inessential especially as energies of these quarks are, on average, lower than energies of quarks evaporated directly. 
Thus, only the neutrino opacity of the photosphere can modify our results. According to a general phylosophy of a 
photosphere's model an opaque "neutrino photosphere" eventually forms p3), but it may happ en too late, when the 



black hole mass is already too small to give noticeable contribution to the integral in Eq. ( |3.11 ) 
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FIG. 7. Constraints on the spectral index n as a function of the reheating temperature Trh from three types of the neutrino 
experiments. 
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A precise calculation of the PBH neutrino background must include also a taking into account of the neutrino 
absorption during a travelling in the space. In the case of the photon background the most important absorption 
process is a pair production on neutral matter [Q due to which the photons from PBHs evaporated earlier than 
z w 700 are absent today. In our case, the analog of an optical depth of the universe for the neutrino emitted at a 
redshift z and having today an energy E is given by the integral 

z 

t(z, E) = cJ a (£7(1 + z')) ■ n{z')^dz' . (5.3) 
o 

Here, cr(E) is the neutrino interaction cross section, n(z) is a number density of the target particles. 

Two processes are potentially " dangerous" : neutrino- nucleon inelastic scattering growing linearly with an energy, 
and annihilations with neutrinos of the relic background: 

v e + N — > e~ + anything, 

(5.4) 

v e + D e (relic) — » >,(/■> + /»)■ (5.5) 

i 

Here, fi are charged fermions (leptons and quarks). As is known, the relic neutrino background exists, with a Planck 
distribution, from the epoch of neutrino decoupling (z « 10 10 ). 

For estimations of the neutrino absorption it is enough to use some characteristic value of a neutrino energy. For 
the neutrino processes in the detectors, discussed in the previous Section, such value is about lOOMeV. In the case 
of the j/A-scattering one has, approximately, 

a vN {E) ~ 10~ 38 • cm- 2 , (5.6) 



= nw „ (io-6 - l(r 7 ) cm- 3 . (5.7) 

(1 + zy 



Using these values and the formula (3.10) for dt/dz, one can easily estimate the integral (5.3) for r(z, E). One obtains 
the important result: 

t vN {z « 10 7 , E = lOOMeF) « 1 . (5.8) 



In the case of the absorption through annihilations of evaporated neutrinos with relic neutrinos the integral ( |5.3| ) 
reduces to |30| ] 

z 

wCI!)* 0.3 -10-" (^).±/(l^)«i>. (5.9) 



This gives 

T ui >(z ~ 2 • 10 6 ,£7 = lOOMeV) < 1 . (5.10) 

From here one can conclude that the neutrino absorption becomes important (for the problems considered in this 
paper) at z ~ 2- 10 6 . It was shown in Sec{[V| that at T RH = 10 10 GeV the redshifts - 10 7 give a maximum contribution 
to the neutrino background (see Fig.ph. B eside s, it was shown that, in general, a redshift value corresponding to such 



a maximum is proportional to T RH (Eq.(iT3)). So, the value z = 2 ■ 10 6 corresponds to Trh ~ 6 • 10 9 GeV , and, 
therefore, for Trh > 6 • 10 9 GeV the neutrino absorption effects are important. 

One should note, at the end, that usually the calculations of these constraints are accompanied by the calculation 
of the bounds based on requirement that the energy density in PBHs does not overdose the universe at any epoch 
(Qbh < 1). For a setting of such bounds one must consider the cosmological evolution of the system PBHs + 
radiation. We intend to carry out these calculations in a separate paper. 
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